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Abstract 

Gauge theories in axial gauges are studied using Exact Renormalisation Group flows. 
We introduce a background field in the infrared regulator, but not in the gauge fixing, 
in contrast to the usual background field gauge. We discuss the absence of spurious 
singularities and the finiteness of the flow. It is shown how heat-kernel methods can 
be used to obtain approximate solutions to the flow and the corresponding Ward 
identities. New expansion schemes are discussed, which are not applicable in covari- 
ant gauges. As an application, we derive the one-loop effective action for covariantly 
constant field strength, and the one-loop /3-function for arbitrary regulator. 
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I. INTRODUCTION 



The perturbative sector of QCD is very well understood due to the weak coupling of glu- 
ons in the ultraviolet (UV) limit, known as asymptotic freedom. In the infrared (IR) region, 
however, the quarks and gluons are confined to hadronic states and the gauge coupling is 
expected to grow large. Thus the IR physics of QCD is only accessible with non-perturbative 
methods. The exact renormalisation group (ERG) provides such a tool UK]. It is based on a 
regularised version of the path integral for QCD, which is solved by successively integrating- 
out momentum modes. 

ERG flows for gauge theories have been formulated in different ways (for a review, see 
0). Within covariant gauges, ERG flows have been studied in fH||6[], while general ax- 
ial gauges have been employed in In these approaches, gauge invariance of physical 
Greens functions is controlled with the help of modified Ward or Slavnov- Taylor identities 
lllJQ'illtni • A different line has been followed in [JTTJ] based on gauge invariant variables, 



e.g. Wilson loops. Applications of these methods to gauge theories include the physics of su- 
perconductors |12| , the computation of instanton- induced effects |l3j , the heavy quark effec- 
tive potential fT3,15|, effective gluon condensation JjiJ, Chern-Simons theory fllTf , monopole 



condensation [18|, chiral gauge theories [jnj, supersymmetric Yang-Mills theories |[2(J, and 



the derivation of the universal two-loop beta function ||21|| . 

In the present paper, we use flow equations to study Yang-Mills theories within a back- 
ground field method. In contrast to the usual background field formalism ||22|| , we use a 
general axial gauge, and not the covariant background field gauge. The background field 
enters only through the regularisation, and not via the gauge fixing. Furthermore, in axial 
gauges no ghost degrees of freedom are present and Gribov copies are absent. Perturba- 
tion theory in axial gauges is plagued by spurious singularities of the propagator due to an 
incomplete gauge fixing, which have to be regularised separately. Within an exact renormal- 
isation group approach, and as a direct consequence of the Wilsonian cutoff, these spurious 
singularities are absent 0. The resulting flow equation can be used for applications even 
beyond the perturbative level. This formalism has been used for a study of the propagator 
for a formulation of Callan-Symanzik flows in axial gauges pi, and for a study of 



Wilson loops 

Here, we continue the analysis of and provide tools for the study of Yang-Mills 
theories within axial gauges. First we detail the discussion of the absence of spurious 
singularities. Then a framework for the evaluation of the path integral for covariantly 
constant fields is discussed. We use an auxiliary background field which allows us to define 
a gauge invariant effective action. The background field is introduced only in the regulator, 
in contrast to the usual background field formalism. This way it is guaranteed that all 
background field dependence vanishes in the infrared limit, where the cutoff is removed. 
We employ heat kernel techniques for the evaluation of the ERG flow. The heat kernel is 
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used solely as a technical device, and not as a regularisation. The flow equation itself is by 
construction infra-red and ultra-violet finite and no further regularisation is required. As 
an explicit application, we compute the full one-loop effective action for non-Abelian gauge 
theories. This includes the universal /^-function at one loop for arbitrary regulator. We also 
discuss new expansions of the flow, which are not applicable for covariant gauges. 

The work is organised as follows. We begin with a brief review of the Wilsonian approach 
for gauge theories. This includes a derivation of the flow equation. We discuss the absence 
of spurious singularities and the finiteness of the flow. This leads to a mild restriction on 
the fall-off behaviour of regulators at large momenta. (Section [Tl]). Next, we consider the 
implications of gauge symmetry. This includes a discussion of the Ward-Takahashi identities, 
the construction of a gauge-invariant effective action, and the study of the background field 
dependence. Explicit examples for background field dependent regulators are also given 
(Section pll) . We derive the propagator for covariantly constant fields, and explain how 
expansions in the fields and heat kernel techniques can be applied in the present framework 
(Section |TV|) . We compute the full one loop effective action using heat kernel techniques. 
We also show in some detail how the universal beta function follows for arbitrary regulator 
functions (Section [V]). We close with a discussion of the main results (Section |VlD and leave 
some more technical details to the Appendices. 



II. WILSONIAN APPROACH FOR GAUGE THEORIES 

In this section we review the basic ingredients and assumptions necessary for the con- 
struction of an exact renormalisation group equation for non-Abelian gauge theories in 
general axial gauges. This part is based on earlier work [7|||. New material is contained in 
the remaining subsections, where we discuss the absence of spurious singularities and the 
finiteness of the flow. 



A. Derivation of the flow 

The starting point for the derivation of an exact renormalisation group equation are 
the classical action Sa for a Yang-Mills theory, an appropriate gauge fixing term 5* g f and a 
regulator term ASk, which introduces an infra-red cut-off scale k (momentum cut-off). This 
leads to a fc-dependent effective action IV Its infinitesimal variation w.r.t. k is described 
by the flow equation, which interpolates between the gauge-fixed classical action and the 
quantum effective action, if ASk and Ff. satisfy certain boundary conditions at the initial 
scale A. The classical action of a non-Abelian gauge theory is given by 

SA[A] = \jd i xF; v {A)F;M) (2-1) 
with the field strength tensor 



2 



F; u {A) = d»K - d u Al + gf\ c A\Al (2.2) 



and the covariant derivative 

D;\A) = 5 a % + gf acb A% [t\ f] = f a bc t a . (2.3) 
A general axial gauge fixing is given by 

S s{ [A} = 1 -Jd^xn^^- 2 n u Al. (2.4) 

The gauge fixing parameter £ has the mass dimension —2 and may as well be operator- 
valued JFJ. The particular examples £ = and £p 2 = —1 are known as the axial and the 
planar gauge, respectively. The axial gauge is a fixed point of the flow @ . 

The scale-dependent regulator term is 

AS k [A, A}= 1 -J d 4 xA; R k f u [A}A b u . (2.5) 

It is quadratic in the gauge field and leads to a modification of the propagator. We have 
introduced a background field A in the regulator function. Both the classical action and 
the gauge fixing depend only on A. The background field serves as an auxiliary field which 
can be interpreted as an index for a family of different regulators -R^A- Its use will become 
clear below. 

The scale dependent Schwinger functional W k [J, A], given by 

expW k [J,A} = JvAexp{-S k [A,A]+ |d 4 a;^j;|, (2.6) 

where 

S k [A, A] = S A [A] + S gi [A\ + AS k [A, A] . (2.7) 
We introduce the scale dependent effective action T k [A, A] as the Legendre transform of 



T k [A, A] = j d A xJ^Al - W k [J, A] - AS k [A, A], A% = 5Wk ^ A] . (2.8) 

For later convenience, we have subtracted AS k from the Legendre transform of W k . Thus 
T k [A, A] is given by the integro-differential equation 

exp -T k [A, A] = jva exp{-^[a] - ^[o] - AS k [a -A,A] + j^ k [A, A\{a-A)}. (2.9) 
The corresponding flow equation for the effective action 
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d t T k [A,A] = ^Ti{G k [A,A}d t R k [A}} (2.10) 

follows from ( |2.9| ) by using (a — A) = 0. The trace sums over all momenta and indices, 
t = Ink. Gk is the full propagator of the field A, whereas A is not propagating. Its inverse 
is given by 

(G tl A,A])-' >,*') = S! ( ^U + R ^ ( 2 ' n » 

There are no ghost terms present in ( |2.10| ) due to the axial gauge fixing. For the regulator 
R k we require the following properties at A — 0. 

fk 2 V 

lim p 2 R k = 0, lim R k ~ p 2 — , (2.12) 

where p 2 is plain momentum squared. Regulators with 7 = 1 have a mass-like infra-red 
behaviour with R k (0) ~ A; 2 . The example in (|3.21| ) has 7 = 1. In turn, regulator with 7 > 1 



diverge for small momenta. The latter condition in (|2.12|) implies that R k introduces an IR 
regularisation into the theory. The first condition in fl2.12|) ensures the UV finiteness of the 
flow in case that G k oc p~ 2 for large p 2 . For covariant gauges this is guaranteed. Within 
axial gauges, additional care is necessary because of the presence of spurious singularities. 
It is seen by inspection of (|2.9|) and ( [2.12| ) that the saddle-point approximation about A 
becomes exact for k — > 00. Here, T k approaches the classical action. For k — > 0, in turn, the 
cut-off term disappears and we end up with the full quantum action. Hence, we confirmed 
that the functional T k indeed interpolates between the gauge-fixed classical and the full 
quantum effective action: 

lim F k [A, A] = S[A] + SJA], (2.13a) 

fc^oo 

lim T k [A, A] = T[A\. (2.13b) 

fe^O 

Notice that both limits are independent of A supporting the interpretation of A as an index 
for a class of flows. It is worth emphasising that both the infrared and ultraviolet finiteness 
of ( |2.10| ) are ensured by the conditions ( |2.12|) on R k . 



B. Absence of spurious singularities 



The flow equation ( |2.10[ ) with a choice for the initial effective action Ta at the initial 



scale A serves upon integration as a definition of the full effective action T = T k=0 . It 
remains to be shown that ( [2.10|) is finite for all k thus leading to a finite T. In particular 



this concerns the spurious singularities present in perturbation theory: the propagator P^ v 
related to S = Sa + Sg is 
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p - irtL 

p2 



n 2 (l + & 2 )p^Pv 



1 {jiyp,, + n v Pv) 



(rip)' 



p* 



np 



(2.14) 



It displays the usual IR poles proportional to 1/p 2 . We observe additional divergences for 
momenta orthogonal to n^. These poles appear explicitly up to second order in l/[np) and 
can even be of higher order for certain (np)-dependent choices of £. For the planar gauge 
^p 2 — — 1, the spurious divergences appear only up to first order. 

This artifact makes the application of perturbative techniques very cumbersome as an 
additional regularisation for these spurious singularities has to be introduced. We argued 
in that these spurious singularities are missing in the flow equation. Here, we further 
the discussion, also providing some information about the intricate limit where the cut-off 
is removed. First of all we derive a bound on the flow (|2.10 ). Then, we argue that this 
bound results in weak constraints on the decay behaviour of the regulator function r for 
large momenta. This is sufficient for providing a well-defined RG flow. 

We start with an analysis of the momentum dependence of the propagator in the presence 
of the regulator. To that end we set the background field to zero, A = 0, and specify the 
regulator as 



COfe 



r(p )p V -r(p )p fl p 



(2.15) 



The IR/UV limits of r, r can be read-off from ( 2.12 ). In ( 2.15|) we did not introduce terms 



with tensor structure (n^p y + n v p^) and n^n u . For the present purpose, the discussion of 



spurious singularities, the choice ( |2.15| ) suffices. Indeed, even r plays no role for the absence 
of spurious singularities in the flow equation approach. The only important term for the 
discussion of spurious singularities is that proportional to the term p 2 S fiy 5 ab . It is this term, 
proportional to the identity operator, that guarantees the suppression of all momentum 
modes for large cut-off. The other tensor structures are proportional to projection opera- 
tors and cannot lead to a suppression of all modes. With a regulator obeying (|2.15| ) the 
propagator takes the form 



D <W , VtiVv . n^p u + n u p^ n^n u 

k.uv — CLl— 5" + «2 1 r «3 57 s r 0-4 , , ? 

pz p^ p z [np) n z p z 



with the dimensionless coefficients 



(2.16) 



a 1 = l/(l + r), 

a 2 = (l + r)(l+^p 2 (l + r))/z. 
a 3 = -(1 + f)s 2 /z , 
a 4 = —(r — f)/z, 



(2.17a) 
(2.17b) 
(2.17c) 
(2.17d) 



and 
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s 2 = (np) 2 /(n 2 p 2 ), (2.17e) 
z = (i + r )[(i + f) s 2 + ( r -f)(i +: p^(i + r ))]. (2.17f) 

Now we evaluate the different limits in p 2 and k important for the approach. To keep things 
simple we restrict ourselves to the case f = and a regulator r leading to a mass-like IR 
behaviour: \im p 2/ k 2^ r(p 2 ) = k 2 /p 2 . For this choice we deduce from ( 2.16Q and fl2.12|) that 
Pk,nv has the limits 

p2/ feoo Pfc - = p - A^o Pk ^ = h + ^ttw 2 ) ' (2 - 18) 



with Pfj, u defined in fl2.14|) . By construction, the propagator (|2.16|) is IR finite for any k > 



Now, the important observation is the following: in contrast to the perturbative propagator 
P^, the limit of P^^v for n P — ► is finite. This holds true even for an arbitrary choice of 
£(p,n) and leads to 

1 + r p 2 (l + r)r p 4 (l+r)(l + ^(l+r))nV l ' J 

Thus is well-behaved and finite for all momenta p. The plain spurious divergences are 

already absent as soon as the infra-red behaviour of the propagator is under control. This 
holds for R with the most general tensor structure as long as it obeys the limits Q2.12|) . For 



example, it is easily extended to non-zero f as long as the regulators r and r have not been 
chosen to be identical. Already in the infrared region f has to be smaller than r in order to 
have a suppression of longitudinal modes at all. So we discard the option of identical r and 
f. 

Still, for np = and large momenta squared y = p 2 the regulator tends to zero and the 
second term in ( 2 . 1 9| ) diverges in the limit y — ► oo proportional to y~ l (r — f) _1 > y d / 2_1 , 



following from ( 2.12Q . Hence, even though the term only diverges for y — > oo, a more careful 



analysis is needed for proving the finiteness of the flow equation. We emphasise that the 
remaining problem is the integration over large momenta in the flow equation and not an 
IR problem at vanishing momentum. Thus, by showing that this problem is absent in the 
flow equation for all k it cannot reappear at k — 0. Indeed, we shall see that finiteness of 
the flow for all k implies a stronger decay of the regulator for large momenta as in ( |2.12| ). 
In turn, one may expect problems for regulators with weaker decay. 



C. Finiteness of the RG flow 

Here, finiteness of the flow equation is proven by deriving an upper bound for the flow 
following a bootstrap approach. The derivation of the flow equation is based on the existence 



of a finite renormalised Schwinger functional for the full theory [21|. In the present context 
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this only implies the existence of a renormalisation procedure for axial gauges, the form 
of which is then determined by integrating the flow. An explicit systematic constructive 
renormalisation procedure is not required. The latter is a problem in perturbative field 
theory: no renormalisation procedure is known, which can be proven to be valid to all 
orders of perturbation theory. 

In the present approach, spurious singularities could spoil the finiteness of Q2.1Q) due to 
infinities arising from the integration of the large momentum domain. For the derivation 
of a bound we can safely assume, that for all k and large momenta p 2 the full propagator 

(2) 

T k is dominated by its classical part (possibly with some multiplicative renormalisation 
constants). Hence for large momenta we can estimate T^\s^ + Sg )~~ l < C[A,A] with 
C[A, A] > 0. Consequently the field independent part of the flow provides a bound on the 
full flow. The only terms that could produce divergences are related to the terms in ( 2.16Q 



proportional to a 2 and a 3 , the source for divergences being z 1 . The coefficient a± of the 
last term in ( |2.16|) also contains z~ l but also an additional factor r. Hence the limit np — > 



can be safely done in the term 04. 

We do not go into the details of the computation. A more detailed derivation and 
discussion is given elsewhere. We quote the result for f = 0. Upon integrating the angular 
s-part of the momentum integration we get an estimate from the part of Tr Pk dtRk with 
the slowest decay for y — > 00 



bound oc 



dyy 



2 VT+^ r\y) 



l + < Jr{y) 



(2.20) 



where the square root terms stem from an integration f_ 1 ds/[s 2 + (1 + £y)r(y)]. Since the 
potential problem only occurs from an integration over large momenta squared y = p 2 , we 
have restricted the y-integral to y > a where a is at our disposal. It can be chosen the 
same for all k. This ensures that the limit k — > can be taken smoothly. The bound ( 2.20Q 



stems from the second term in (|2.16|) proportional to a 2 . Eq. (|2.20|) is finite for regulators 
r that decay faster than y~ 5 . Without spurious singularities, r has to decay stronger than 
y~ 2 , see (|2 . 1 2| ) . Hence we have a mild additional constraint due to the fact that the full 
propagator G fc does not introduce an additional suppression. Typically, the regulator is 
chosen to decay exponentially for large momenta. Similar finite integrals as in Q2.20| ) also 
occur in field dependent terms in the flow, as we shall see later in Sect. M. 

This analysis shows the finiteness of the flow ( 2.10|) and supports the claim that the flow 



equation provides a consistent quantisation procedure for gauge theories in axial gauges. The 
bound also marks the use of Callan-Symanzik (CS) type flows (Rk oc k 2 and r(y) oc y^ 1 ) 
as questionable in axial gauges. It has been already mentioned in H] that such a choice 
requires an additional renormalisation. The presence of contributions from all momenta at 
every flow step makes the limit k — > an extremely subtle one. This limit is very sensitive 
to a proper fine-tuning. In axial gauges, this problem for CS flows gets even worse due 
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to the spurious singularities. We know that a consistent renormalisation procedure in the 
axial gauge is certainly non-trivial. For CS-type flows, one is back to the original problem 
of spurious singularities in perturbation theory, but with a more difficult propagator and 
additional renormalisation problems. A recent calculation of perturbative corrections to the 
Wilson loop has indeed shown that formulations in axial gauge with a mass term for the 



gauge field meet problems p5| , p6| . The massless limit of this observable did not coincide 
with the well-known result. In turn, for regulators which decay faster than r(y) ~ y~ 5 , the 
problem is cured. 



III. SYMMETRIES 

In this section, the issue of gauge invariance of physical Greens functions, controlled by 
modified Ward-Takahashi identities, is studied. We discuss the role of background fields, 



which, in contrast to the usual background field method will only be introduced for 

the Wilsonian regulator term. The Ward-Takahashi identities for the quantum and the 
background field are derived. We define a gauge-invariant effective action as it follows from 
the present formalism, and discuss its background field dependence. Finally, we discuss the 
background field dependent regularisation. 

A. Modified Ward-Takahashi Identities 

We now address the issue of gauge invariance for physical Greens functions. The problem 
to face is that the presence of a regulator term quadratic in the gauge fields is, a priori, in 
conflict with the requirements of a (non-linear) gauge symmetry. This question has been 
addressed earlier for Wilsonian flows within covariant gauges [|],|)]|||7||| . The resolution 
to the problem is that modified Ward-Takahashi identities (as opposed to the usual ones) 
control the flow such that physical Greens functions, obtained from at k = 0, satisfy the 
usual Ward-Takahashi identities. 

The same line of reasoning applies in the present case even though in the presence of the 
background field A some refinement is required || . In this particular point it is quite similar 



to the symmetry properties of the full background field formalism as discussed in |10[ . The 
background field makes it necessary to deal with two kinds of modified Ward-Takahashi 
Identities. The first one is related to the requirement of gauge invariance for physical Green 
functions, and is known as modified Ward Identity (mWI). The second one has to do with 
the presence of a background field A in the regulator term Rk, and will be denoted as the 
background field Ward-Takahashi Identity (bWI). 

To simplify the following expressions let us introduce the abbreviation 8^ and 8^ for the 
generator of gauge transformations on the fields A and A respectively: 
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8^A = D(A)u 
8„A = 



S u A = 

8 U A = D(A)u. 



(3.1a) 
(3.1b) 



The action of the gauge transformations 8^ and 8^ on the effective action T k can be com- 
puted straightforwardly. It is convenient to define 



W h [A, A; u] = 6 u T h [A, A] - Tr (n^u) —n v A v + -Tr uo [G k [A, A], R k [A] 
W k [A,A; u] = S u T k [A,A] - ^Tru [G k [A, A], R k [A] 



(3.2a) 
(3.2b) 



In terms of ( |3.2| ), the behaviour of r^LA, A] under the transformations 8^ and 8^, respec- 
tively, is given by 








Wk[AA;u] 
m[A,A;uj} 

Eq. ( |3.3b| ) is valid for regulators R k that transform as tensors under 8^, 

S u R k [A] = \R k [A],u . 



(3.3a) 
(3.3b) 



(3.4) 



Eq. ( |3.3a| ) is referred to as the modified Ward-Takahashi identity, and ( |3.3b| ) as the back- 
ground field Ward-Takahashi identity. 

Let us show that (|3.3|) is consistent with the basic flow equation ( 2.10|) . With consistency, 
we mean the following. Assume, that a functional T k is given at some scale k which is a 
solution to both the mWI and the bWI. We then perform a small integration step from k 
to k' — k — Ak, using the flow equation, and ask whether the functional T k i again fulfils the 
required Ward identities ( |3.3| ). That this is indeed the case is encoded in the following flow 
equations for (|3.3|), namely 



d t W t [A,A;uj] 
d t W k [A,A;u,] 



2^ \ k 
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H^-w g ^-a 



W k [A,A;u] 



(3.5a) 
(3.5b) 



■ a ji-t. Eq. (B.5I) states that the flow of mWI is zero if the 
mWI is satisfied for the initial scale. The required consistency follows from the fact that the 
flow is proportional to the mWI itself ( |3.5a|) , which guarantees that (|3.3a|) is a fixed point 
of ( |3.5a|) . The same follows for the bWI by using ( |3.5b|) . There is no fine-tuning involved 
in lifting a solution to ( |3.3a ) to a solution to (|3.3b|) . It also straightforwardly follows from 
( gja| ) and 03351) . 

We close with a brief comment on the use of mass term regulators. Such a regulator 
corresponds simply to R k = k 2 and leads to a Callan-Symanzik flow. The regulator is 
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momentum-independent, which implies that the loop term in ( |3.2a| ) vanishes identically. 
Hence one concludes that the modified Ward identity reduces to the usual one for all scales 
k. This happens only for an axial gauge fixing J/J. 



B. Gauge invariant effective action 

Returning to our main line of reasoning and taking advantage of the results obtained 
in the previous section, we define a gauge invariant effective action only dependent on A 
by identifying A = A. It is obtained for a particular choice of the background field, and 
provides the starting point for our formalism. 



It is a straightforward consequence of the mWI ( fj.3a|) and the bWI ( |3.3b| ) that the effec- 



tive action Tk[A, A] is gauge invariant - up to the gauge fixing term - under the combined 
transformation 

(<L + S u )T k [A,A] = Tr n^d^-^-rivA,. (3.6) 

We define the effective action T k [A] as 

t k [A]=T k [A,A = A]. (3.7) 
The action T k [A] is gauge invariant up to the gauge fixing term, to wit 



S u t k [A] = Tr { n ^u)^-n v A v \ . (3.8) 



This follows from ( |3.6| ). Because of ( [2.13b|) , the effective action f fc =0 [A] is the full effective 



action. The flow equation for T k [A] can be read off from the basic flow equation ( 2.10|) 



d t t k [A] = ~Tr {G k [A, A] d t R k [A}} , (3.9) 

Notice that the right-hand side of ( |3.9|) is not a functional of I^LA]. The flow depends on 
the full propagator G k [A, A], which is the propagator of A in the background of A taken at 
A = A. Thus for the flow of T k [A] one needs to know the flow (of a subset) of vertices of 
5 2 r k [A, A]/(5A) 2 at A = A. Still, approximations, where this difference is neglected are of 
some interest [PTfl . 

We argue that ( |3.8| ) has far reaching consequences for the renormalisation procedure of 
T k [A] as is well-known for axial gauges and the background field formalism. T k [A] is gauge 
invariant up to the breaking due to the gauge fixing term. We define its gauge invariant 
part as 

T k , im [A] = T k [A] - S g f[A] (3.10a) 
S u rk,h*[A] = . (3.10b) 
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Eq. ( |3.10| ) implies that the combination gA is invariant under renormalisation, dt(gA) = 0. 
If one considers wave function renormalisation and coupling constant renormalisation for A 
and g respectively 

A -> Z l J 2 A (3.11a) 
9 - Z a g (3.11b) 

we conclude that 

Z g = Z F 112 . (3.12) 



C. Background field dependence 

By construction, the effective action r^LA, A] at some finite scale k ^ will depend on 
the background field A. This dependence disappears for k = 0. The effective action fk[A] is 
the simpler object to deal with as it is gauge invariant and only depends on one field. As we 
have already mentioned below (|3.9| ), its flow depends on the the propagator <5^r fc [A, A] at 
A = A. Eventually we are interested in approximations where we substitute this propagator 
by 5\T k . The validity of such an approximation has to be controlled by an equation for the 
background field dependence of T k [A, A]. The flow of the background field dependence of 
T k [A, A] can be derived in two ways. Sjd t T k can be derived from the flow equation fl2.10| ), 



-Ld t Y k [A,A] = i^Tr {G k [A,A}d t R k [A}} . (3.13) 

The flow dt5^T k follows the observation that the only background field dependence of 1^ 
originates in the regulator. Thus, 5^T k is derived along the same lines as the flow itself and 
we get 

d t -tv k [A,A] = ±Tr d t ^G k [A,A]^R k [A}X, (3.14) 

which turns out to be important also for the derivation of the universal one loop /3-function 
in Sect. V B . The difference of ( |3.13| ) and ( 3.14 ) has to vanish 



[^,dt]T k [A,A) = 0. (3.15) 

Eq. ( |3.15| ) combines the flow of the intrinsic A-dependence of 1^ [A, A] (|3.14|) with the A- 
dependence of the flow equation itself ( |3.13| ). It provides a check for the validity of a given 
approximation. Using the right hand sides of ( |3.13| ) and ( |3.14| ) the consistency condition 
( P-15| ) can be turned into 
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Tr 



sr (2) 

Gk ^ Gk d t Rk 



Tr G k "-jf G k d t T k 



(2) 



(3.16) 



where 



T k 2) [A,Ar(x,x') 



ab 



S 2 T k [AA] 
5A%(x)5Al{x>) ' 



(3.17) 



With ( |3.16|) , we control the approximation 



S 2 T k [A,A] 



5 A 5 A 



A=A 



5 2 T k [A] 
5 A 5 A 



+ sub-leading terms 



(3.18) 



For this approximation the flow (|3.9|) is closed and can be calculated without the knowledge 



of T k 2 \ but with f k 2 \ Amongst others, the approximation ( 3.18Q is implicitly made within 



proper-time flows, where the use of heat-kernel methods is even more natural ||28|| . This is 
discussed in p9j (see also [27|). Let us finally comment on the domain of validity for the 
approximation ( |3.18|) . In the infrared k — > 0, the dependence of the effective action I^LA, A] 
on the background field A becomes irrelevant, because the regulator i?fc[A] tends to zero. 
Therefore we can expect that ( |3.18| ) is reliable in the infrared, which is the region of interest. 



D. Regulators 

We have seen that the symmetries of the effective action T k and the flow crucially depend 
on the properties of i?&[.A], in particular the construction of a gauge invariant effective action. 



The regulator has to transform as a tensor under gauge transformations of A, (3.4). Here we 



specify a general class of regulators which has this property and is well-suited for practical 



applications. As already argued in section IIB, the infrared regularisation is provided by r 



whereas r ^ only gives different weights to the longitudinal degrees of freedom, see ( 2.1 5|) . 
In the following we set f = 0. We choose 

R k [A]=D T r{D T ) (3.19) 

with the yet unspecified function r. We introduced Dt, the Laplace operator for spin 1, 

D^M ■= -{D P D p ) a \A)5, u - 2gF;l{A) (3.20) 

and Dt = Dt(A). For vanishing background field the Laplacean Dt reduces to the free 
Laplacean -Dt(O) = p 2 . In this case we have R k = p 2 r(p 2 ). Written in terms of some general 
Laplace operator P 2 (A), a typical example for the regulator functions R k (P 2 ) and r(P 2 ) is 

^( p2 ) = — pfL i > r ( p2 ) = — J/i.2 i ( 3 - 21 ) 

expF 2 /k' s — 1 expI JZ /k z — l 
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which meets the general properties as described in ( |2.12| ). Eq. ( |3.21| ) is an example for a 
regulator with a mass-like IR behaviour, 7 = 1. More generally the IR/UV conditions for 
Rk in fl2.12p translates into 

2 




lim U- r = 0, lim r ~ — (3.22) 




for the function r. 



IV. ANALYTIC METHODS 

In this section we develop analytical methods to study flow equations for gauge theories 
in general axial gauges. The flow equation is a one-loop equation which makes it possible 
to use heat kernel techniques for its solution. The main obstacles, technically speaking, 
are the constraint imposed by the modified Ward identity and the necessity to come up 
with a closed form for the full propagator. We first derive such an expression for the case 
of covariantly constant fields within general axial gauges. In addition a generic expansion 
procedure in powers of the fields is discussed. Finally, we give the basic heat kernels to be 
employed in the next section. 



A. Propagator for covariantly constant fields 

We derive an explicit expression for the full propagator for specific field configura- 
tions. This is a prerequisite for the evaluation of the flow equation fl2.10|) . To that end 
we restrict ourselves to field configurations with covariantly constant field strength (see 
e.g. 01), namely D^F up = 0. This is a common procedure within the algebraic heat 
kernel approach. We also use the existence of the additional Lorentz vector to demand 
n^A^ = n^F^ = 0. That this can be achieved is proven by the explicit example of n p = S^o 
and (Ap) = (A = 0,Ai(x)). These constraints lead to 

[D li ,F up ]=0, (4.1a) 
n p A p = (4.1b) 
n^ u = 0. (4.1c) 

To keep finiteness of the action of such configurations we have to go to a theory on a finite 
volume. However, the volume dependence will drop out in the final expressions and we 
smoothly can take the limit of infinite volume. For the configurations satisfying (|4.1| ) we 
derive the following properties 

[D\D p ] = -2gF, p D p , (4.2a) 
D T , W D P = —D^D 2 , (4.2b) 
[n p D p , D p ] = 0. (4.2c) 
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Defining the projectors P n and Pd with 

Pn.nu = — 7T~ , (4.3a) 
tit 

Pd iM , = A^A, (4.3b) 

we establish that 

P D D T = -P D D 2 P D , P n D T = ~P n D 2 (4.4) 

holds true. After these preliminary considerations we consider the gauge-fixed classical 
action given in ( |2.1|) .We need the propagator on tree level to obtain the traces on one-loop 
level. The initial action reads 

T A [A] = S A + S gi . (4.5) 
From (|4.5| ) we derive the full inverse propagator as 



rg?[4 A] = ( D^ u + (D fl D u ) ab + ^~ 2 n,n u S ab ) + 0(g 2 ; D T , D,D V ) . (4.6) 



The inverse propagator ( fO| ) is an operator in the adjoint representation of the gauge group. 
We now turn to the computation of the propagator ( 2.11 ) for covariantly constant fields. 
Using O, (0) and §3), we find 

n \& A-\ ah ( ( ai \ n a2 n , a 3 n . n % , n tl a 4l n v \ , , 

Gk[A A] ^ = ~ (fej^ + D ^ D " + n ^15HnD) + D »m^ + ) ' (4 - 7) 

with the dimensionless coefficient functions 

1 

ox = — , (4.8a) 

l ~T} l V D) (s 2 + rn[l-D>ai + r D ) ] y 1 , (4.8b) 
(l + r D ) ^ l 

(s 2 + m[l-D 2 Z(l + r D )]) , (4.8c) 



a 2 



{1 + rn, 



(s 2 + r J ,[l-D 2 e(l + r D )]) 1 . (4.8d) 



Notice that a\ is a function of while 02, 03 and 04 are functions of both D 2 and (nD) 2 . 
We also introduced the convenient short-hands 

r T = r k (D T ), r D = r k (-D 2 ) } s 2 = (4.9) 

The regulator function, as introduced in ( |3.19| ), depends on D T . The dependence on D 2 , as 
apparent in the terms a 2 , a 3 and a 4 , comes into game due to the conditions ( |4.1| ) and (fOl). 
They imply 
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r k (D T )D„ = D^r k (-D 2 ), r k (D T )n„ = n^-D 2 ) , (4.10) 

which can be shown term by term for a Taylor expansion of r k about vanishing argument. 
For vanishing field A = the propagator (Op reduces to the one already discussed in M. 
There, it has been shown that the regularised propagator ( [4.7] ) (for r ^ 0) is not plagued by 
the spurious propagator singularities as encountered within standard perturbation theory, 
and in the absence of a regulator term (r = 0). For the axial gauge limit £ = the expression 
( |4.7|) simplifies considerably. With ( |4.6| ) and ( |4.9| ) we get 

r \A] -( 1 ^ -d - p i n » nD n 

^[A\ { DT{1+rT) )^ "*D*(l + r D )(* + r D r v+ n*D*{l + r I ,)(*' + r D ) L ' t ' 

nD n v ro 

+ "D^l + r D )(s 2 + r D ) + D 2 (l + r D )(s 2 + r D ) ' (4 ' 11} 



The propagators (4/7) and (|4.11j) are at the basis for the following computations. Notice that 



this analysis straightforwardly extends to approximations for r^LA, A] beyond the one-loop 

,[IV [ 



level. Indeed, it applies for any Tk[A, A] such that T k '[A, A] is of the form 



inD £nD 



Here, the scale-dependent functions f^ T and fJ? D can depend on Dt, D 2 and nD. In 
turn, the functions fj? D and fj} n can depend only on D 2 and nD. An explicit analytical 
expression for the full propagator, similar to (|4.7|) , follows from (|4.12|) . Such approximations 



take the full (covariant) momentum dependence of the propagator into account. The inverse 
propagator ( |4.6|) corresponds to the particular case f DT = f DD = 1, f nD = 0, and f nn = 



B. Expansion in the fields 



Even for analytic calculations one wishes to include more than covariantly constant 
gauge fields, and to expand in powers of the fields, or to make a derivative expansion. 
Eventually one has to employ numerical methods where it is inevitable to make some sort 
of approximation. Therefore it is of importance to have a formulation of the flow equation 
which allows for simple and systematic expansions. 

In this section we are arguing in favour for a different splitting of the propagator which 
makes it simple to employ any sort of approximation one may think of. For this purpose 
we employ the regulator Rk[D 2 (A)]. This is an appropriate choice since it has no negative 
eigenvalues. We split the inverse propagator into 

r?j?[4 = A£ - (2gF£ - (D,D u r b ) (4.13) 
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with 



A J = {-£ 2 (1 + r D )} ab 5, u + ^n^. 



(4.14) 



The operator A can be explicitly inverted for any field configuration (and A — A). We have 

1,1 1 



A" 1 



-11 + 



D 2 (l + r D ) D 2 (l + r D )l + £D 2 (l + r D 
With ( f4.13p and ( |4.15| ) we can expand the propagator as 

oo 

G k [A, A] = A- 1 [(2gF -D®D) A" 1 



-Pr,.. 



(4.15) 



(4.16) 



n=0 



where (D ® D) a ^ u {x,y) = D a ^D^5(x — y). For £ = (the axial gauge), A 1 can be neatly 
written as 



A (£ = 0) 



1 



£> 2 (l + r D ) 



;n--Pn), 



(4.17) 



which simplifies the expansion ( |4.16| ). The most important points in (|4.16 ) concern the 
fact that it is valid for arbitrary gauge field configurations and each term is convergent for 
arbitrary gauge fixing parameter £. Moreover such an expansion is not possible in the case 
of covariant gauges. Both facts mentioned above are spoiled in this case. 



C. Heat kernels 



We present closed formulae for the heat-kernel of the closely related operators Dt and 
— D 2 = Dt + 2gF. These are needed in order to evaluate the traces in ( 5.15|) . We define 
the heat-kernels as Kq{t) = exp{rO}(x,x) 



K D 2(T) 



K. 



(27T) 4 

2rF K D ,(r 



(4.18a) 
(4.18b) 



where X M = ip^ + D^ in the corresponding representation. Here we used that 2gF commutes 
with X M for covariantly constant fields. All kernels are tensors in the Lie algebra ( K^d t is 
also a Lorentz tensor because of the prefactor). For the calculation of the momentum integral 
we just refer the reader to the literature (e.g. p0| ) and quote the result for covariantly 
constant field strength 



K D 2{t) 



167T 2 T 2 



det 



rgF 



sinh rgF 



1/2 



K_ Dt (t) =exp(2rgF) K D 2 (r) 



(4.19a) 
(4.19b) 
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Here, the determinant is performed only with respect to the Lorentz indices. For the com- 
putation of the one-loop beta function we need to know K(t) in ( 4. 19| ) up to order F 2 
(equivalently to order r°). Expanding K D 2 in rgF we get 

Km(r) = -L- (1 - ^g 2 (F\ p ) + 0[r, (gF) 3 }. (4.20) 



16tt 2 Vr 2 12 y v 7PP 

With ( ^20|) and the expansion (exp2rgF)^ = 1 + 2rgF flv + 2r 2 g 2 (F 2 ) flu + 0[r, (#F) 3 ] we 
read off the coefficient of the K(r) proportional to F 2 , 

Tr Kdi I pi = - J-j jWsUfA] , (4.21a) 
1 2D 

Tr F_ Dt | f2 = — - 2 -Ng 2 S A [A] , (4.21b) 

where the trace Tr denotes a sum over momenta and indices. We have also used that 
Sa[A] — | J" tr f F 2 with tr ft a t b = —\5 ab . Since the operators Dt and D 2 carry the adjoint 
representation the trace Tr includes tr a ^ with 2Ntift a t b = tr a d£ a i b . 

V. APPLICATIONS 

In order to put the methods to work we consider in this section the full one-loop effective 
action for SU(N) Yang-Mills theory which entails the universal one-loop beta function for 
arbitrary regulator function. 

A. Effective action 

For the right hand side of the flow we need 

T k [A,A]= l -Jz F {t)ti f F 2 {A) + S gf [A}+0[{gAf ) g 2 dA} ) tr f t a t b = -U ab (5.1) 

where tr^ denotes the trace in the representation R, R = f stands for the fundamental rep- 
resentation, R = ad for the adjoint representation. Only the classical action can contribute 
to the flow, as n-loop terms in (|5.1|) lead to n + 1-loop terms in the flow, when inserted on 
the right hand side of (3J). This Ansatz leads to the propagator ( 4.1 1|) which together with 



our choice for the regulator ( |3.19| ) is the input in the flow equation (|3.9|). We also use the 



following in the evaluation of the different terms in 

tr D 2 = 4tr D ® D (5.2) 

With this we finally arrive at 

f _ 1 I d t r(D T ) 1 d t r(-D 2 ) 1 d t r(-D 2 ) \ 
1 k ~ 2 \ 1 + r{D T ) 2 1 + r{-D 2 ) + 4 s 2 + r{-D 2 ) J ' { > 
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where the trace Tr contains also the Lorentz trace and the adjoint trace ti a d in the Lie 
algebra. The first term on the right-hand side in ( |5.3| ) has a non-trivial Lorentz structure, 
while the two last terms are proportional to 5^ u . We notice that the flow equation ( |5.3| ) 
is well-defined in both the IR and the UV region. We apply the heat-kernel results of 
section [IV Q to the calculation of (|5~3f) . To that end we take advantage of the following 



fact: Given the existence (convergence, no poles) of the Taylor expansion of a function f(x) 
about x = we can use the representation 

f(-O) = f(-d T )ex P {rO}\ T=0 (5.4) 

Due to the infrared regulator the terms in the flow equation ( |5.3| ) have this property, where 
O = Dt, D 2 . Hence we can rewrite the arguments Dt and — D 2 in ( |5T3| ) as derivatives 
w.r.t. r of the corresponding heat kernels K^o T ij) and K D 2{t). Applying this to the flow 
equation (|5.3| ) we arrive at 

d t r{-d T ) 1 d t r(-d T ) 

■ TtK- Dt (r) - -— — , Q , TrA D2 (r) 



d t t k = i 



l + r (-d T ) ' 21 + r(-d T 

(5.5) 



+ - fd Pn (P 2 n- d r) d tr(p 2 n -d T ) ^ 2 ^ R 

4 J n p 2 n + (pi - d T )r(p 2 n - d T ) ° 



J T = 



The two terms in the first line follow from ( |5.3| ). The last term is more involved because 
it depends on both D 2 and nD due to s 2 = (nD) 2 /n 2 D 2 . We note that nD = (nd) holds 
for configurations satisfying ( |4.1a ) and only depends on the momentum parallel to n 



A'' 



Furthermore it is independent of the gauge field. Now we use the splitting of (p M ) = (p n ,p) 
where p n = P n p and p = (1 — P n )p. The heat kernel related to D 2 follows from the one for 
D 2 via the relation Kg 2 (r) = z ^K D 2(t) as can be verified by a simple Gaufiian integral in 
the p n -direction. 

With these prerequisites at hand, we turn to the full effective action at the scale k, which 
is given by 

where A is some large initial UV scale. We start with the classical action = Sa + Sg. 
Performing the /c-integral in (|5.6|) we finally arrive at 

T k [A] = (l + S (f - 7(1 - 7)) Ink/ A) S A [A] 

oo 

+Sgf[A] + C m (k 2 /k 2 ) AT^[gF/k 2 } + const. (5.7) 



m=l 



The combination Sa + Sgf on the right-hand side of ( |5.7| ) is the initial effective action. 
All further terms stem from the expansion of the heat kernels ( |4.19[ ) in powers of r. The 



terms ~ r 2 give field-independent contributions, while those ~ r 1 are proportional to 
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trF and vanish. The third term on the right-hand side of ( |5.7| ) stems from the r° coef- 
ficient of the heat kernel. This term also depends on the regulator function through the 
coefficient 7 ( |3.22| ). All higher order terms ~ r m ,m > are proportional to the terms 
C m (k 2 /A 2 ) AT^[gF/k 2 ]. These terms have the following structure: They consists of a 
prefactor 



C m (x) 



1 



4m ml 



x 



(5.8a) 



and scheme-dependent functions of the field strength, A T^ m '[gF], each of which is of the 
order 2 + m in the field strength gF. They are given explicitly as 



AT^[gF] = B% TrK%{0) + (sf + B£>) TrK$(p) . 



(5.8b) 



Here, i^ } (0) and i^ T (°) denote the 

expansion coefficients of the heat kernels. We use 
the following identity 



'(in) 



/M(o) = f(d T y 



lr=0 ' 



(5.9) 



and f^ m '(x) = (d x ) m f(x). In addition, the terms in (|5.8b ) contain the scheme-dependent 
coefficients 



m 

Dfl 2 

B1 D 



i + n 



--b Dt 
2 m 



(m) 



(0) 



— dx [d x - -ad a 



ri(x) 



ri(x)Jr x {x) + a 



(5.10a) 
(5.10b) 

(5.10c) 



The coefficients B Dt , B d and B nD follow from the first, second and third term in 
We introduced dimensionless variables by defining r\{x) = r(xk 2 ) and fx(x) = dtri(x) = 
—2xk 2 r'(xk 2 ) = —2xr[(x), in order to simplify the expressions and to explicitly extract the 
^-dependence into ( p.8a| ) . The explicit derivation of B nD is tedious but straightforward and 
is given - together with some identities useful for the evaluation of the integral and the 
derivatives - in appendix |A|. All coefficients B Dt , B d2 and B nD are finite. The appearance 
of roots in the coefficient B nD is not surprising after the discussion of the absence of spurious 



singularities in section II B 



In particular, we can read off the coefficients for m = which add up to the prefactor 
of the classical action in ( |5.7| ): 



Bn 



27, 



-7, 



T>nD 

Bq 



7), 



(5.11) 



where we have used ( |A.5| ) in the appendix. Together with the heat kernel terms proportional 
to r° given in QC2ID this leads to (|577|). 
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This application can be extended to include non-perturbative truncations. The flow of 
the coefficients (|5.8b| ) becomes non-trivial, and regulator-dependent due to the regulator- 
dependence of the coefficients ( [5.101 ). Then, optimisation conditions for the flow can be 
employed to improve the truncation at hand [31] . 

Finally, we discuss the result (|5.7f) in the light of the derivative expansion. Typically, 
the operators generated along the flow have the structure F fk[{F> 2 + k 2 )/A 2 ] F, and similar 
to higher order in the field strength. For dimensional reasons, the coefficient function fk(x) 
of the operator quadratic in F develops a logarithm ~ lnx in the infrared region. An 
additional expansion of this term in powers of momenta leads to the spurious logarithmic 
infrared singularity as seen in (|5.7[ ). To higher order in the field strength, the coefficient 
function behave as powers of 1/(D 2 + k 2 ), which also, at vanishing momenta, develop a 
spurious singularity in the IR, and for the very same reason. All these problems are absent 
for any finite external gluon momenta, and are an artifact of the derivative expansion. 
A second comment concerns the close similarity of ( |5.7|) with one-loop expressions found 
within the heat-kernel regular isat ion. In the latter cases, results are given as functions of the 
proper-time parameter r and a remaining integration over dlnr. Expanding the integrand 
in powers of the field strength and performing the final integration leads to a structure as 
in ( |5.7| ), after identifying r ~ k~ 2 . In particular, these results have the same IR structure 
as found in the present analysis. 



B. Running coupling 

We now turn to the computation of the beta function at one loop. We prove that 
the result is independent of the choice of the regulator and agrees with the standard one. 
However, it turns out that the actual computation depends strongly on the precise small- 
momentum behaviour of the regulator, which makes a detailed discussion necessary. 

Naively we would read-off the /^-function from the t-running of the term proportional to 
the classical action Sa in (|5.7|). Using Q3.12| ) leads to d t \nZ g = — ^dtlaZp. We get from 



(E3 

**=(f- 7(1-7))^* "> d t \nZ g = -(f -1(1-1)) ^- 2 +0(g% (5.12) 



We would like to identify (3 = dtlnZ g . This relation, however, is based on the assumption 
that at one loop one can trade the IR scaling encoded in the t-dependence of this term 
directly to a renormalisation group scaling. This assumption is based on the observation that 
the coefficient of Sa[A] is dimensionless and at one loop there is no implicit scale dependence. 
It is the latter assumption which in general is not valid. A more detailed analysis of this fact 



is given in [21]. Here, we observe that the background field dependence of the cut-off term 



inflicts contributions to d t Z F S c \. These terms would be regulator-dependent constants for 
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a standard regulator without A. As mentioned below (|2.5|), one should see the background 
field as an index for a family of different regulators. We write the effective action as 

r k [A, A] = r M L4] + r k , 2 [A] + r kt3 [A, A] . (5.13) 

The second term only depends on A and is solely related to the A-dependence of the regula- 
tor. The last term accounts for gauge invariance of T k under the combined transformation 
S w + 6 W . This term vanishes in the present approximation, because of the observation that 
our Ansatz is invariant - up to the gauge fixing term - under both 5^ and 5^ separately. 
The physical running of the coupling is contained in the flow of r^ifA]. This leads to 

(3 = -\d t Z F + ±d t Z F}2 , (5.14) 

where Z F>2 is the scale dependence of T k >2 oc Z f ^Sa[A]. We rush to add that this procedure 
is only necessary because we are interested in extracting the universal one-loop /3-function 
from the flow equation. For integrating the flow itself this is not necessary since for k = 
the background field dependence disappears anyway. For calculating d t In Z F2 we use ( |4.11| ) 
and ( |5.2j ) and get 



5A° 1 2 l \D T + R k [D T }5A- 2 —D 2 + R k [—D 2 ] 5A^ 

R' k [-D 2 ] 6D 2 \ 



(5.15) 



A(-nD) 2 + R k [-D 2 ] <L4£j' 
where we have introduced the abbreviation 

R' k {x) = d x R k {x). (5.16) 



For the derivation of (|5.15 ) one uses the cyclycity of the trace and the relations ( f4.2| ). We 



notice that Q5.15 ) is well-defined in both the IR and the UV region. The explicit calculation 



is done in appendix 0. Collecting the results (|B.2j ),( [Blf ),( [B.4|) we get 



d t 5 A T k [A, A = A]\ F 2 = -^f- 2 7(1 - 7 ) 6 A S A [A] - d t Z F , 2 = 7(1 - 7 ) (5.17) 

l07T^ \bll Z 



We insert the results ( p,12j ) for d t Z F and ( |5.17|) for d t Z F ^ 2 in ( p. 14j ) and conclude 



/3 = -^^4 + 0(/). (5.18) 

which is the well-known one- loop result. For regulators with a mass-like infrared limit, 7 = 1, 
there is no implicit scale dependence at one loop. It is also worth emphasising an important 
difference to Lorentz-type gauges within the background field approach. In the present case 
only the physical degrees of freedom scale implicitly with t — Ink for 7 / 0. This can be 
deduced from the prefactor 7(1 — 7) in ( p . 1 7| ) . Within the Lorentz-type background gauge, 
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this coefficient is -g-(l — 7) ||2T|. The difference has to do with the fact that in the axial 
gauge one has no auxiliary fields but only the physical degrees of freedom. In a general 
gauge, this picture only holds true after integrating-out the ghosts. This integration leads 
to non-local terms. They are mirrored here in the non-local third term on the right hand 
side of the flow (|5.5|) and in the third term on the right hand side of (|5.15|) [see also (B.4)]. 



VI. CONCLUSIONS 

We have shown how the exact renormalisation group can be used for gauge theories 
in general axial gauges. We have addressed various conceptual points, in particular the 
absence of spurious singularities and gauge invariance, which are at the basis for a reliable 
application of this approach. We have shown that spurious singularities are absent provided 
that the regulator decays stronger than (p 2 )~ 4 for large momenta. In turn, regulators 
with milder decay are highly questionable. At least they are subject to a renormalisation 
of the flow itself, which implicitly brings back the problem of spurious singularities. This 
concerns in particular the mass regulator = k 2 , see also [||. 

Our main goal was to develop methods which allow controlled and systematic analytical 
considerations. The formalism has the advantage that ghost fields are not required. Also, 
no additional regularisation - in spite of the axial gauge fixing - is needed. This is a 
positive side effect of the Wilsonian regulator term. In addition, we worked in a background 
field formulation, which is helpful in order to construct a gauge invariant effective action. 
Also, it allows to expand the flow equation around relevant field configurations. Instead of 
relying on the standard background field gauge, we have introduced the background field 
only in the regulator term. The axial gauge fixing is independent on the background field. 
This way, it is guaranteed that the background field dependence vanishes in the IR limit. 
It is important to discuss how this differs from the usual background field approach to 
Wilsonian flows. In both cases, applications of the flow require an approximation, where 
derivatives w.r.t. the background field are neglected, cf. ( |3.18| ). In the present approach, this 
approximation improves in the infrared, finally becoming exact for k = as the background 
field dependence disappears. For the background field gauge this does not happen, because 
the full effective action still depends non-trivially on the background field. 

As an application, the full one-loop effective action and the universal beta-function have 
been computed. This enabled us to address some of the more subtle issues of the formalism 
like the implicit scale dependence introduced by the cutoff, which has properly to be taken 
into account for the computation of universal quantities, and the scheme independence of 
the beta-function. The equation which controls the additional background field dependence 
introduced by the cutoff contains the related information. 

These results are an important step towards more sophisticated applications, both nu- 
merically and analytically. A natural extension concerns dynamical fermions. The present 
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formalism is also well-adapted for QCD at finite temperature T, where the heat-bath singles- 
out a particular Lorentz vector. Here, an interesting application concerns the thermal pres- 
sure of QCD. 
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A. EVALUATION OF THE ONE LOOP EFFECTIVE ACTION 



The calculation of the last term in ( |5.7|) is a bit more involved. Note that the following 
argument is valid for m > — 1, m > — 1 is of importance for the evaluation of ( |5.7| ), m = — 1 
will be used in Appendix [Bl We first convert the factor r m+1 / 2 appearing in the expansion 
of the heat kernel using T l l 2+m = (— l) m+1 -2= J dzd™ +1 e~ Tz2 . We further conclude that 



m D - — [j n J~ (Pn-d T )d t r(p 2 n -d T ) +1 2 



:-i) 



m+1 



471 

'_]\m+\ 



dp n dz d. 



(m+1 



d t r(p 2 n - d T 



471 



dp n dzd\ 



m+1 



Vl + (Pi ~ d T )r(p 2 n - d T 
d t r(z 2 +p 2 n ) 



r=0 



2 „2 

Vn 

z 2 +pl 



+ r(z 2 + p 2 



(A.l) 



The expression in ( |A.1[ ) can be conveniently rewritten as 
1) 
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nD 



\m+l ^oo ^2tt 

- / dx I i 

8n Jo Jo 

.^m+l ^oo ^ [ 



d. r 



\ m+1 

-ad n 
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Jo 



dx d T ad n 



x 



x 

m+1 



d t r(x) 



a sin <p + r ( x ) 
d t r(x) 



a=l 



r(x)Jr(x) + a 



(A.2) 



a=l 



where x = z 2 + p 2 n and sin 2 (ft = Pn/(z 2 + pi). It is simple to see that — (l/x)ad a is a 
representation of d z i on sin 2 (f> = p 2 J (z 2 +p 2 ) and d x a representation of d z i on functions of 
x only. The expression in ( |A.2j ) is finite for all m > 0. Evidently it falls of for x — > oo. For 
the behaviour at x = the following identity is helpful: 



i \ m+1 m+1 

d x --ad Q ^ ^ \ ni + 1— i 

x 



EC- 1 ) 7 



i=0 



(A.3) 



Eq. ( A.3 ) guarantees that the integrand in ( |A.2j ) only contains terms of the form 
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a H7F7TT7 (l + rr) '~""7 (A ' 4) 

with i — 0, m+ 1. For a; — > one has to use that <9 t r — > 2nr and r — > — . The terms of 
integrand in ( |A.2j ) as displayed in ( |A.4j ) are finite for x — 0. 

We are particularly interested in Bq D relevant for the coefficient of Sa in the one loop 
effective action ( pTTQ . With ( [/Cfl it follows 

1 r j_ f a 1 .* A ^(x) 



/| / 

4 Jo 



£™ D = __ / dx d r - -ad, 



x ' \Jr(x)\jr(x) + a 



a=l 




- 2 4£L] -4(1-7), (A.5) 



r(x)y / l + r(x) \A + 



x=0 



where we have used dtr(z) = —2zd z r(z) and the limits for dtr(z — > 0) = 272; 7 ,r(^ — > 0) 
00) = 0. 



B. ^-DERIVATIVES 



For the calculation of ( |5.15 ) the following identity is useful: 



Tr ( ^r°) eT ° = - T ^ eT ° ^ (B.i) 

where we need (|B.1|) for O = D 2 and O = —D?- Now we proceed in calculating the first 
term in ( |5.15| ) by using a similar line of arguments as in the calculation of Q5.7D and in 
Appendix |A|. We make use of the representation of r _1 = J °° dzexp —rz and arrive at 

\rd ( R ' k[DT] Wt ) - X Trd ( 1 6 K M \ 

2 Tldt \D T + Rk [D T }jA^) ~ 2^ \-d T + R k [-d T ) t5A» K - D ^ T) ) t = 

1 r dx f R' k [x] \ Ng*20 6 

^-^fd-7)^(^] + OM). (B.2) 

Note that dt acts as —2xd x on functions which solely depend on x/k 2 . The term R'/ (1 + r) 
is such a function. The second term can be calculated in the same way leading to 

1 a f —R' k [D 2 } 5 -1 1 f»dx a ( R' k [x] \Ng 2 4 6 



f 9 *(l^p^]^ )-li T 9 ' (.iTTRJ 16^3M« 
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The calculation of the last term in (|5.15| ) is a bit more involved, but boils down to the same 
structure as for the other terms. Along the lines of Appendix |A| it follows that this term 
can be written as 

'Trd f 5D2 \ -\rdl fdv gMl^d ^ 6 K lW ] 

8 Tldt \(-nD)i + R k[ -Di}SA«j ~ 8 Tl % [J d *V B + R k \vl ~ d T ] ^TSA*^ > 



T=0 



^5( i -^4 (s ^ 1+0[sl) (B - 4) 

Note that when rewriting the left hand side of ( |B.4| ) as a total derivative w.r.t. A this also 
includes a term which stems from j^{nD) 2 . This, however, vanishes because it is odd in p n . 
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